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Pade-Borel approximation of the continuum limit of strong coupling lattice fields: 
Two dimensional non- linear 0{N) sigma model at > 3 
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Based on the strong coupling expansion, we reinvestigate two dimensional 0{N) sigma model by the use of 
Pade-Borel approximants. The conventional strong coupling expansion of the mass square M in momentum 
space in j8 = l/g~ is inverted to give j8 expanded in l/M. Borel transform of j5 with respect to M is carried 
out and the result is improved as the rational function by Pade method. We find the behavior of Pade-Borel 
transformed bare coupling at 18th order is consistent for N >3 with that of continuum scaling to the four-loop 
perturbation theory. We estimate non-perturbative mass gap at A' > 3 and find the agreement with the exact 
result by Hasenfratz et.al. 
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I. INTRODUCTION 

Nearly four decades ago, the quark confinement was shown 
by Wilson at the strong bare coupling region Ol. For weak 
coupling, perturbation theory clarified for the Yang-Mills sys- 
tem that bare coupling g tends to vanish as the lattice spacing 
fl ^ 101 ■ The motivation of the present work is to attempt to 
extrapolate the large a behavior of bare coupling to the asymp- 
totically free behavior at weak coupling. For the purpose, we 
like to reformulate the strong coupling expansion by changing 
the primary variable from bare coupling to the lattice spacing 
itself. 

Lattice serves us a suitable regularization, since in lattice 
field theories the lattice spacing a explicitly appears in the 
action and enters into the physical quantities. For instance, 
the dimensionless correlation length ^ represents a physical 
length scale divided by a. It is given at strong coupling as a 
series ^ in (g^)^^ and it determines, in an implicit manner, 
the a dependence of the bare coupling. Mutual roles of g and 
^ are exchanged by inverting the relation Thus we ad- 
dress the question whether the small £, series of g allows us 
to confirm directly the weak coupling behavior predicted by 
perturbation theory. 

In the present paper, in non-linear 0{N) model at two di- 
mensions, we make an attempt to approximate the asymptotic 
behavior of bare coupling in the continuum limit via its large a 
expansion. The model is of interest as a testing ground of our 
approach, since it enjoyes asymptotic freedom and dynamical 
mass generation for > 3 13[|. In addition to the large limit, 
we also consider the case of finite A^. 

As the basic variable, rather than the correlation length in 
lattice space, we adopt mass M in momentum space defined 
by the zero momentum limit of the two-point field correla- 
tion. The lattice spacing a is included in the mass which is 
rescaled to be dimensionless and then M must vanish in the 
continuum limit a ^ (see ( |2.7| i). Now the strong coupling 



expansion gives series of M in j3 = {g^)^^- By inverting the 
series, we express )3 as a power series in 1 /M, which is equiv- 
alent to large a expansion. As it would be, naive series fails 
to confirm the continuum behavior of j5 . However, it is non- 
trivial and interesting to examine, when both Pade and Borel 
techniques are applied on the series, whether the continuum 
scaling emerges at finite A^ or not. 

Before proceeding to following sections, we remark the role 
of Borel transform in our approach. We use Borel transform 
as a device of dilation operation around the continuum limit. 
The response of scale transformation on f{M) is probed by 
rescaling M to AM in / and taking the X ^ limit. Then, it 
is said / scales with the exponent A if 



/(AM) ^ A '^/(M) , A 0. 



(1.1) 



The above criterion of scaling is implemented by introducing 
8 defined by 



A = l-5, 0<5<1 



(1.2) 



and performing expansion to some finite orders in 5 10. Ht] . 
Suppose the function approaches to M^^. Then, expand- 
ing it to 5^ and setting 5 = 1, one has (M(l — 5))^^ — > 
rTTTxSi-ATuM^^- Further, if we take the limit M °°,L ^ 
°o with M/L = M fixed, we obtain 
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That is, the limit A — > has transitioned to the limit 5 — > 1 
with the cut off L. Then the scaling behavior with exponent 
A manifests itself in the power of M(= M/L). Note the uni- 
versal quantity A is left unchanged. On the other hand, when 
same operation is acted on / in the series form ^fl/tM^'^ valid 
at large M, we have / = Y,{ai^/kl)M^'' with a larger conver- 
gence radius, which is just the Borel transform of the original 
series. We thus interpret the Borel transform as a realization 
of scale transformation. We do not need integrating / back to 
/. Though the information of /(M) over the whole range of 
M is not obtained, what we need in lattice field theories is the 
behavior of /(M) in the neighborhood of M = 0. 
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II. DESCRIPTION OF THE MODEL 

On the two-dimensional square lattice, the continuous spin 
fields d = (cJi , C52, • • • , On) are set on every sites. The action 
of the system is given by 



n 11 = 12 

where ei = (1,0), 62 = (0, 1) and 



i3 = 



1 



(2.1) 



(2.2) 



The fields are constrained to satisfy at every sites, a = A^. 

The mass variable M defined via the zero momentum limit 
of the propagator (Ln^xpl'P ' ' ^''(ii))) ' is given by 



2DX 



(2.3) 



where susceptibility x and second moment fi are, respectively, 
given by ^Lni^m ■ ^(n)) and ^ = ]^I„n2((j(0) ■ 
cj(n)). D denotes the dimension of lattice space and D ~ 2 
in the present work. Let us summarize the continuum limit of 
the model and large a expansion of j5 . 

The perturbative renormalization group predicts that, for 
N >3, the correlation length behaves at weak coupling as 



^ = exp 



N-2 



1/{N- 



L Qk 

k=\ 



(2.4) 



where the multiplied constant Q is specified only non- 
perturbatively. Hasenfratz et. al. has computed it via ther- 
modynamic Bethe ansatz |@], giving 



,1-71/2. _L 



N-2 



(2.5) 



The terms flAjS^'' (k ^ 1,2,3, ■■ ■) in (12.4b are contributions 
of k + 2-loop levels and three- |0] and four-loop fsl] results 
were computed in the literature. They are given as 



1 



02 



NiN - 2) 
1 



(-0.0490- 0.014 W), 



A?2(A?-2)2 

X (0.0444 + 0.0216A^+0.0045A^2_o.0129A^^). (2.6) 

Though three and higher loop contributions disappear in the 
continuum limit for the bare coupling, we cannot take out the 
limit because only the series to finite order is at hand. Hence, 
we include known three- and four-loop contributions in our 
analysis. 

It is known thatM^'/^ has functional form of j3, the same as 
(12.4b but with another multiplicative constant, say Cm- How- 
ever, Monte Carlo data showed that the difference is less 
than a percent at = 4. Since the two constants agree with 
each other in the large limit, the difference between and 



Cm niay actually be negligible for all > 3. Thus the estima- 
tion of the mass gap via strong coupling expansion becomes 
the estimation of C^ and this is one of the aims of our work. . 

Since the mass M approaches to ^ tjjg continuum Umit, 
the physical mass of dimension 2 is given by 



'■phys 



(2.7) 



where Ai is the finite mass scale given by 
-27rA^j3i f2nN^\-ViN-2) 



Al = fl 2 exp 



A^-2 



\N-2J 



k=l 



(2.8) 

From ( |2.4b and ^ M^'/^^ have continuum /3 to four-loop 
order as a function of M, 



N-2 



4nN 



1 /I X 



-27rA?(A?-2)fli+log(ilog^ 



kN{N -2)\og-i, 



nN{N-2Y{\ogPf 
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where 



1 



(2.9) 



(2.10) 



On the series expansion at large M, we borrow the result 
in the work of Butera and Comi ifToll who computed strong 
coupling series of x and /i to jS^'. Using the result, we have 
expansion of M in powers of j3 , 



M 



1 



2(3+2A0„ 2(16 + 32A^+17A^2_^2A^3) ^3 



2+N 



{2+Ny-{4+N) 



{2+Nf 

By inverting the above relation, we have 

2(10A^+19y 8(14A^ + 25)x4 



(2.11) 



A^ + 2 A^ + 2 

2 (338A'3 + 2593A^2 ^ 6084A^ + 45 12) 
{N + 2f{N + A) 



(2.12) 



Based upon the series ( 12.12b . we discuss the approximation 
of the continuum limit by the use of Pade-Borel approxima- 
tion scheme. We attempt to recover the asymptotically free 
behavior ( 12. 9b from ( 12.121 ) and then estimate C^ . 
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III. LARGE LIMIT 

Large limit serves us a good bench mark of our approach. 
So we consider that case first and then turn to finite in the 
next section. 

In the large limit, only the one-loop contribution to j3 
survives to give 



1 

47t 



log(x/C^), ^ (32)-^/2 ^ 0.17677669 



(3.1) 

As briefly presented in the introduction, Borel transform is 
given by a certain limit of delta expansion |5i]. Explicitly, 
the logarithm is expanded and gives at 5 = 1 that log(ji:/(l — 
5) ) — ?> log jc + Yd= 1 T ^'^ order L. Then using the asymptotic 
expansion ^4=1 V = logi + 7£ +0(L^') {Je denotes Euler's 
constant), we have logx log(jcL) + Ye in the L ^ °° limit. 
Let X be small enough with x = xL kept finite, then the result 
represents Borel transform of logx. Denoting the operation 
of Borel transform by ^ we thus find ,'^[logx] = logx + Ye- 
Using abbreviated simbol j3 = .^[j3], we then obtain 



l5^1-(log{x/Cl) + YE)^Pa 



The large M expansion of j3 reads 

j3 = X - + 20^-3 _ 112/ + 676x^ - 4304x^ - 



(3.2) 



(3.3) 



Borel transform of the above series results to divide the nth 
order coefficient by the factorial of «, 



2!" 



20 
3!"' 



x^- 



112 

j 

4! 



676 



x'- 



4304 

i 

6! 



(3.4) 

Then as a crucial step, we use Pade method to extrapolate the 
above series to larger x. The resultant Pade-Borel approxi- 
mants enable us to capture the scaling behavior to be seen in 
the scaling region as we can see below. 

As a preliminary study, we have examined the behaviors of 
[m/n] approximants of j3 over almost possible pairs of m,n at 
orders m + « = 4, 5 , ■ • • , 20. On the contrary to the condensed 
matter models undergoing second order phase transition, crit- 
ical behavior of the present model is known from perturbation 
theory as logarithmic and slowly varying. Hence it is con- 
ceivable that good behaviors come from the cases where the 
difference between m and n is small. The numerical experi- 
ment confirmed it is indeed the case. We have also compared 
the approximants of three types, Pade-Borel, Borel only and 
Pade only improvements. The result at 6th order is shown in 
FIG. 1. As already reported in jsl, Borel only improvement 
is not sufficient for observing the asymptotic freedom. Pade 
only case (jSp/sj in FIG. 1) is also insufficient as is clear from 
FIG. 1 . However, Pade-Borel approximant shows enough im- 
provement for quantitative approximation. Though Pade only 
approximation is found to be improved at higher orders, the 
best performance is achieved by Pade-Borel approximant at 
every order we analyzed. We therefore focus on Pade-Borel 
approximant hereafter. 




— P[3/3] 
~ P[3/3] 

~ P[6/0] (Borel) 

- - 1-Ioop (Naive) 
---- 1-Ioop (Borel) 



FIG. 1 : Plot of improved j3 and jS at 6th order. Two broken lines (one 
for j8 and the other for j8) represent behaviors at continuum. Hori- 
zontal axis corresponds to logx = log(l/M) and logx = Iog(l/M) 
(for Pade only case). 



Now, let us turn to the evaluation of the mass gap by esti- 
mating . Since we know information at weak coupling, the 

estimation is carried out by fitting Pcont to [m/n] order approx- 
imants of j3, j3[,„/„], by adjusting the value of Q. In practice, 

we consider the difference between j3[„,/„] and Pamt and plot 
the difference by changing the value of . At just proper 
value of Q , the two functions touch with each other at a point 
xo and the difference is tiny over an interval including xq. A 
typical case is shown in FIG. 2 and the result of estimation of 
Q is shown in Table |T] For the reason previously written, we 
list only the results around the diagonal Pade. Though the 




logx 

FIG. 2: Subtracted function jSp/g] - Pcont = j8[9/9] - {4^(logx/C| + 
Ye)}- Plotted curve is for = 0.17868 

reason is not known to us, the orders 6, 10, 14 and 18 give the 
best approximation among nearby orders. 



IV. FINITE DOWN TON = 3 

In this section we study the weak coupling behavior from 
Pade-Borel approximants for a finite number of spin compo- 
nents. First we discuss Borel transform of ( 12.91 ) to compare it 
with Pade-Borel approximants of large M series ( |2.12| i. 

Let us consider Borel transform of the two-loop contribu- 
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TABLE I: Evaluation result denoted as Capp of the non-perturbative 
constant in ^(log(x/C|) 
(32)-'/- =0.17677669- ■•. 



Je)- Rigorous value of 



[m/ri\ 


^app 


[m/n] 


Capp 


[3/3] 


0.18327 


[7/7] 


0.17911 


[4/3] 


0.18734 


[8/7] 


0.18041 


[3/4] 


0.18722 


[7/8] 


0.18038 


[4/4] 


0.18463 


[8/8] 


0.17972 


[5/5] 


0.18138 


[9/9] 


0.17868 


[6/5] 


0.18266 


[10/9] 


0.17900 


[5/6] 


0.18264 


[9/10] 


0.17901 


[6/6] 


0.18178 


[10/10] 


0.17875 



tion. We find 

[log log x] = log log X 



Ye , C(2) -rl 



logx 2(logx)^ 



0((logx)-3). (4.1) 



The result is simplified by absorbing Je into the log. Then we 
obtain 



?[loglogx] -log(logx + 7£)- 



C(2) 



2(logi + 7£)2 



-O((logx-) 



(4.2) 

Note that the second term should be included when the four 
loop contribution is taken into account. At two- and three- 
loop orders, we need only the first term. In a similar manner, 
for contributions at three- and four-loop orders, we find as 
follows: 



1 



1 



logx 
logx 

(l0gx)2 

loglogx 
(loglogx)^ 

^ (logx)2 J 



logx + 7£ 
log(logx + 7£) 
logx+7£ 
1 

(logx + 7£)2 

log(logx + 7£) 
(logx + 7£)2 

(log(logx + 7g))2 
(logi+7£)2 



0((logi)-3), 
+ 0((logx-)-3), 
0((logx)-3), 
+ 0((logx-)-3), 



(4.3) 
(4.4) 
(4.5) 
(4.6) 



-0((logx)-3). (4.7) 



Thus the result of Borel transform to the four-loop level 
reads, 

A^ — 2 X 1 1 X 
d = (log-T + 7£)H log[-(log^ + 7£)] 

-2nN{N~2)ai+\og[{{log^ + YE)] 

+ J 

nN{N-2){log^ + rE) 



TtN{N-2f{\og^ + YEf 
AnNiN - 2)(-fli +Nnal - 2nNa2) 



-2(l+27rA^(A^-2)fli)log[l(log^ 



-Ye)] 



1 X 

-iog[2(iog^ + r£)]^ 



C(2) 



4;rA^(log-^ + 7£)^ 









1-loop 

2-loop 






3-loop 

4-Ioop 


-4 -2 2 4 6 






-4 -2 2 4 G 


-4 -2 


2 4 G 



N=20 



N=8 



N=6 



N=5 



N=4 



N=3 



(4.8) 



FIG. 3: Plots of jSp ygj at A' = 20, 8,6,5,4,3 and ficont at one and two- 
loop results (plus three and four-loop results fox N = 3,4 and 5) as 
functions of logx. 
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At large M, we have from (12.121 1. 



2(10A^+19)x^ 8{14N + 25)x 



3!(A^ + 2) 



4!(A^ + 2) 



(4.9) 

As in the previous section, we further improve the large M 
series by Pade method. We have checked that also at finite 
A^, diagonal Pade provides best behaviors. Skipping low order 
results, we explicitly present only the results at 18th order for 
various A^. FIG. 3 shows the plots of jSp /gj and ^com at one- 
and two-loop levels (at = 3,4,5, ^cont at three- and four- 
loop levels are also plotted) as functions of logx. At > 6, 
three- and four-loop ^com are very close to that at two-loop at 
X > and we have omitted them. At A^ = 3,4, though the scal- 
ing to the four-loop level is not so clear, the behavior of jSp /gj 
is roughly consistent with the continuum one for logi G [2,4]. 
At A^ = 5, linear-like behavior with correct slope is observed 
around logx ~ 2, which signals scaling behavior From A^ ~ 8, 
we observe continuum scaling at the two-loop level. 

Now, having examined continuum scaling, we evaluate 
constant as in the same manner at A^ = oo. Namely we 

consider jSp /9] — ^cont and search for the value of Q by fitting 

^cont to j3[9 /9] by changing values of . The result is sum- 



TABLE II: Result of estimation of the constant (implied by Capp)- 
The last column shows the result of Botera and Comi llOr. 



A' 


^2— loop 
^app 


^3— loop 
^app 


^—loop 
^app 




Cbc 


3 


0.0068 


0.0094 


0.0112 


0.0125 




4 


0.0336 


0.0373 


0.0398 


0.0416 


0.039 


5 


0.0584 


0.0615 


0.0639 


0.0652 


0.065 


6 


0.0771 


0.0797 


0.0818 


0.0826 


0.084 


7 


0.0913 


0.0934 


0.0953 


0.0955 




8 


0.1021 


0.1038 


0.1055 


0.1054 


0.106 


9 


0.1106 


0.1121 


0.1136 


0.1132 




10 


0.1175 


0.1187 


0.1201 


0.1195 


0.121 


11 


0.1231 


0.1242 


0.1255 


0.1247 




12 


0.1278 


0.1288 


0.1299 


0.1290 


0.130 


13 


0.1318 


0.1326 


0.1337 


0.1327 




14 


0.1352 


0.1359 


0.1369 


0.1358 


0.137 


15 


0.1381 


0.1388 


0.1397 


0.1386 





marized in TABLE II and FIG. 4. We may say that, for all 
N >3, especially for N ^ 3 and 4, j3[9/9] yields a good four- 
loop estimation of non-perturbatice constant . However we 
see that, at A^ > 8, the estimated value Capp is slightly larger 
than the exact one. Note that an excess of the estimation is 
observed also in the large A^ limit. On the contrary, for N <7, 
Capp is slightly smaller than C^. We like to discuss the is- 
sue in the next section. To summarize, we conclude that the 
approximation level is satisfactory. 



. 16 
. 14 
. 12 



. 08 
. OS 
. 04 
. 02 



FIG. 4: Plot of true value of (solid curve) and its estimation at 
A' = 3,4,5,6, 15 (black points) carried out via Pade-Borel ap- 
proximants jSp/gj compared with four-loop ficont- 



V. DISCUSSION 

From previous two sections, at larger A^, we found that the 
four-loop estimation of gives excess to the exact value. For 
example, at [9/9] approximants, the excess reads ~ 0.0001, 
- 0.001 1 and - 0.0019 at A^ 8, 15 and oc, respectively. As 
long as X is not so small, every multiple-loop (above one-loop 
level) contribution decreases as A^ becomes large and vanishes 
at A^ = oo. Hence, at large enough A^, the five, sixth, — loop 
contributions may be safely neglected in our study. Then the 
main factor of the discrepancy would come from the lattice 
artifact. To find evidence, let us discuss the large A^ limit since 
that case provides us quantitative example as we can see be- 
low. 

At N = °°, j3 behaves for small M as 



1 M 

log 

47t ^32 



M , M 

(log — 

327r^ ^32 



l) + 0(M2logM). (5.1) 



Note that this is just the one-loop result. The second and 
higher order terms represent lattice artifacts which disappear 
in the continuum limit. They involve the logarithm and de- 
lay the approach of j3 to the continuum limit. In fact, use of 
Borel transform has the notable advantage that it reduces the 
correction to 



'^ = -4^^^°^32-^^^-m 



Here we have used 



^[MlogM] = -M, 
SS\M\ = 0. 



(5.2) 



(5.3) 
(5.4) 



For original j3 (see (15. Il l), the second term is of orderMlogM, 
but for /3, M and the deviation from the asymptotic scaling 
is much reduced when M is small enough. The correction, 
however, still affects the small M behavior of the transformed 
bare coupling. We have examined scaling and evaluated 

by keeping the first order correction From TABLEllIll 

it is apparent that incorporation of 0{M\ogM) = 0{a^ log a^) 
term improves the approximation. At [9/9] order the excess 
is only ^ 0.0001. This means that Pade-Borel approximants 
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actually recover the small M behavior very well, but in the 
same time, the residual effect of the correction is still non- 
negligible for higher accuracy. We thus find that main factor 
of the discrepancy comes from the lattice artifact as long as 
is large enough. 



TABLE IIL Result of estimation of the constant (implied by Capp) 

in the large A' limit when the correction — to the asymptotic scal- 
ing is taken into account. Only the results of diagonal approximants 
are shown. Rigorous value of is (32) "'/^ = 0.17677669 - -. 



[m/n] 


Capp 


[2/2] 


0.187126 


[3/3] 


0.177873 


[4/4] 


0.178509 


[5/5] 


0.177355 


[6/6] 


0.177484 


[7/7] 


0.176892 


[8/8] 


0.177036 


[9/9] 


0.176891 


[10/10] 


0.176900 



Next, consider the case of lower N = 3 ^ 1 where the es- 
timated value of is slightly lower than the exact one. As 
a typical case, consider the N = 3 case. Three- and four-loop 
effects contribute to pcom ^ 0.4694 at x = xq (log.Jo ~ 3.2) by 
amounts ^ 0.02 and ^ 0.01, respectively. Though they carry 
with small fractions of total j3, they are not negligible at all, 
since Ccpp increases by 0.0026 and 0.0018 when three- and 
four-loop effects is taken into account, and the magnitude of 
itself is small. Therefore, loop contributions above four 
would be still active for estimating and even have possibil- 
ity to push Capp be larger than . For small A^, in addition to 
the lattice artifact, a discrepancy may come also from lack of 
higher loops. 

On the lattice artifact, it is crucial to reduce it for obtaining 
precise result for all A^ > 3. It has reported in fill] that the 
standard action gives a^(logfl^)^/('^^^' as the leading lattice 
artifact near the continuum limit. It has the maximum value 
at A^ = 3 giving contribution ^ a^(logfl^)^. Borel transform 
would reduce the effect of such a logarithmic term but the 
effect would remain to obscure the asymptotic scaling at finite 

X. 

In general, the leading lattice artifact may not be known 
completely. Then, one way to resolve the issue is to construct 
or use lattice action in which such artifacts are reduced from 
the outset. As an example, we report the result of Symanzik's 
modification of lattice action fl2ll in the large A^ limit. In 
Symanzik improvement program, one generalizes the action 
element from 1 - (7n ■ C7n+e^ to Lf=o^A:5'n • ^n+ken ■ By expand- 
ing the action in a and minimizing lattice artifact at the level of 
action, one can obtain the optimized set of coefficients {k = 
0, 1,2,3, • • • ,K). Then, the direct effect is the modification of 
the unperturbed propagator from [2^^=12(1 — cosap^)]^' = 
E|ifl^/'^ "L|ifl'*P^i/12 + 0(aV^)]"' to the one closer to the 



continuum limit [a^Y.^p^] For instance, to the first order 
(K ~ 2) we have 

5 8 1 
~ 3 '^osflp^ + - cos2flp^)]"' 

= [fl25^/_^6^^6/90 + O(oV)]-'. (5.5) 

At the infinite order {K = °o) the action becomes infinite series 
composed of field couplings between two sites along jj. of all 
distances. The result in momentum space is simple modifi- 

TABLE IV: Ratio of Capp (approximant of C^ ) to the exact value of 
in standard, first order and infinite order improved actions. The 

blanks represent absence of extremum zero of P^„f„] — jS . However, 
even in those cases, the difference (the subtracted function) exhibits 
almost stationary behavior around the point, say also xq, at which 
P[n/n] ~ P vanishes and the first derivative takes minimum value. The 
estimation of C^ at such xq yields accurate values. 



[m /n] 


standard 


1st order 


infinite order 


[3/3] 


1.03673 


1.01102 




[4/4] 


1.04443 


1.01137 


0.999244 


[5/5] 


1.02604 


1.00505 




[6/6] 


1.02830 


1.00516 


0.999914 


[7/7] 


1.01320 


1.00168 




[8/8] 


1.01665 


1.00199 


0.999990 


[9/9] 


1.01077 


1.00103 




[10/10] 


1.01145 


1.00103 


0.999999 



cation of the propagator to the continuum limit [fl^L/j f ^] 
Since in the large A^ limit, j3 is given by the gap equation writ- 
ten only with the propagator with mass square M, we can eas- 
ily obtain the large M series both at first- and infinite-order 
improved actions (For detailed presentation, see the first ref- 
erence in 101). For example, at the first order it follows 



1 



-n (2;r)2 
1 5 



1157 8419 



rCOS/?^+ gCOs2/9^) 



-0{M-^) 



36M^ 36M4 



(5.6) 



At infinite-order improvement, the right-hand side becomes 
just the integral of (A^ + Ln=i,2/'^) ' expansion of j3 
in 1/M is straightforward. It now suffices for us to repeat 
the same procedure for the approximation of ( 13.1b and the 
constant at the first and inifinite orders of improved ac- 
tions. Here note that the change of action induces the change 
of the value of non-universal C^. = 0.2377607 - ■■ and 
0.2851456-- - at first and infinite orders, respectively. TA- 
BLE |IV] summarizes the result of our approximation. The 
improved action improves the approximation accuracy both 
at the first and at infinite orders. Though the improved lattice 
action is conventionally used in the Monte Carlo analysis and 
perturbation theory, it is also useful in our approach. 
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In the present work, we have analyzed Pade-Borel approx- 
imants of strong coupHng expansion in non-Hnear a model 
and have found good behaviors approximating the continuum 
limit. We close the paper by pointing out that, even work- 
ing with the standard action, further higher order computa- 
tion would improve the result for all including the limit 



N ^ °°. Pade-Borel approximants may become effective at 
larger x (smaller M) and the two unwanted effects, lattice arti- 
facts and omitted loop contributions, would be weaker there. 
Then, continuum scaling at smaller M with a clearer sign of 
asymptotic freedom near N ^3 would be seen, which allows 
us accurate evaluation of the mass gap for all > 3. 
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